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Abstract. We calculate the width for the lu — > rnr decay in nuclear matter in a hadronic model including 
mesons, nucleons and A isobars. We find a substantial width of the longitudinally polarized lu modes, 
reaching ~ lOOMeV for mesons moving suitably fast with respect to the nuclear medium. 

PACS. 25.75.Dw relativistic heavy-ion collisions - 21.65.+f nuclear matter - 14.40.-n mesons 



The dilepton measurements in the CERES JjJ and HE- 
LIOS Q experiments have indicated that the masses and/or 
widths of light vector mesons undergo large modifications 
in nuclear matter. Clearly, since the mesons interact strongly 
with the medium, this fact is not at all surprising. In- 
deed, in-medium modifications of hadron properties are 
predicted in a variety of theoretical calculations [H,f§,B|Ji. 
||,|,0,0,0[l|,|l|,|l|,|l6) (for a recent review see |[L7, 
lqj). An interesting factor brought in by the presence 
of the medium is that processes which are forbidden in 
the vacuum by symmetry principles are now made pos- 
sible. The constraints of Lorentz-invariance, G-parity, or 
isospin invariance in isospin-asymmetric media pS|,pO[, 
are no longer effective. An example of such an "exotic" 
phenomenon which becomes possible and significant in 
the presence of nuclear matter is the decay of u) — ► rnr. 
This process is, apart for small isospin-violation effect, 
forbidden in the vacuum. Q In this paper we show that 
the matter-induced width for this process is large. For lu 
moving with respect to the medium with a momentum 
above ~ 200MeV (such momenta are easily accessible in 
heavy-ion collisions) the corresponding width, at the nu- 
clear saturation density, is of the order of lOOMeV. In ad- 
dition, we find very different behavior of the longitudinally 
and transversely polarized lu mesons, with the former ones 
being much wider than the latter ones. 

Our calculation is made in the framework of an effec- 
tive hadronic theory. Mesons interact with the nucleons 
and A isobars, and the interactions are assumed to have 



the usual form used in many other calculations and fits. 
We work to the leading order in the nuclear density. Based 
on the experience of other in-medium calculations we hope 
that this approximation should be sufficient up to densities 
of the order of the nuclear saturation density. To this lead- 
ing order only the diagrams shown in Fig. 1 contribute. In 
these diagrams the nucleon lines include the propagation 
of occupied states of the Fermi sea. The "bubble" dia- 
gram (a) has been analyzed by Wolf, Friman, and Soyeur 
in Ref. where the role of the lu — a mixing mecha- 
nism has been pointed out. In this process the lu meson is 
first converted, via interaction with the nucleons, into the 
scalar-isoscalar a meson, which in turn decays into two 
pions. The relevance of "triangle" diagrams (b) has been 
shown out in Ref. Note that in any formal counting 
scheme (low density, chiral limit, large number of colors) 
the diagrams (a) and (b) are of the same order and consis- 
tency requires to include both. Our present calculation of 
lu — > 7T7T in nuclear matter includes a further contribution 
of diagrams (c-d) with the Z\(1232) isobar. Among other 
possible resonances, the A is the most important one due 
to the large value of the 7r7VZ\ vertex and small A — N 
mass splitting. 

The solid line in Fig. 1 denotes the in-medium nucleon 
propagator, which can be conveniently decomposed in the 
free and density parts |23| 



G(k) = G F (k) + G D (k) 
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In the vacuum the partial width for the decay uj — » tt + it~ 
is only ~ 0.2MeV, and is due to the small isospin breaking and 
the resulting p — lu mixing. In this paper we are not concerned 
with this negligible effect. 



where k is the nucleon four-momentum, M denotes the nu- 
cleon mass, Ek = V M 2 + k 2 , and k F is the Fermi momen- 
tum. The diagram (a) is non-zero only when one propaga- 
tor is Grj, and the other one G F . The only non- vanishing 
contributions in diagram (b) involve one Gd propagator 
and two G F propagators. Diagram (a) involves the inter- 
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We use g u — 9. For the ttNN vertex we use the pseu- 
doscalar coupling, with the coupling constant g^NN = 
12.7. The same value is used for ga-NN- The airir cou- 
pling constant is taken to be equal to g a7TT! = 12.8 m^, 
where m n = 139.6MeV is the physical pion mass (this 
value follows from the fit done to tttt scattering phase 
shifts done in Ref. ^l|). The ttNA vertex has the form 
VnNA = (/irJV A/fn^jP^T- ', where p^ is the pion momen- 
tum, T is the g — ► I isospin transition matrix, and the 
coupling constant J^na — 2.1 ptjPl 

The amplitude, evaluated according to the diagrams 
depicted in Fig. 1 (a-d) can be uniquely decomposed in 
the following Lorentz-invariant way: 



M = ^(App + Bu^ + Cg M ), 



(5) 



Fig. 1. Diagrams contributing to the uj — > tttt amplitude in 
nuclear medium. The incoming u has momentum q and po- 
larization e. The outgoing pions have momenta p and q — p. 
Diagrams (b-d) have corresponding crossed diagrams, not dis- 
played. The Feynman rules are given in the text. 



mediate cr-meson propagator, which we take in the form 
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Here the mass and the width of the a meson are chosen in 
such a way that they reproduce effectively the experimen- 



tal tttt scattering length at q 



(780MeV) 2 , which 



is the relevant kinematic point for the process at hand. 
From this fit we find m a = 789MeV and r a = 237MeV. 
Note that m w and m a are very close to each other, which 



enhances the amplitude obtained from diagram (a) 21 



The double line in diagrams (c-d) denotes the A prop- 
agator 
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This formula corresponds to the usual Rarita-Schwinger 
definition p4|,p5| with the denominator modified in order 
to account for the finite width of the A resonance, r A — 
120MeV. 

We assume that the ujNN and uiAA vertices have the 
form which follows from the minimum-substitution pre- 
scription and vector-meson dominance applied to the nu- 
cleoli and the Rarita-Schwinger p3] Lagrangians: 



V 

v 



NN 
IJ.a/3 



AA 



9u [~ 



(4) 



Possible anomalous couplings can be incorporated at the 
expense of having more parameters. The results presented 
below do not depend qualitatively on the form of the cou- 
pling, as long as it remains strong. The coupling constant 
g u can be estimated from the vector dominance model. 



where p is the four- momentum of one of the pions, q 
is the four-momentum of the uj meson, u is the four- 
velocity of nuclear matter, and e specifies the polarization 
of uj. Our calculation is performed in the rest frame of nu- 
clear matter, where u = (1, 0, 0, 0). In this reference frame 
the amplitude j\4 vanishes for vanishing 3-momentum q, 
as requested by rotational invariance. Hence, the process 
uj — > tttt occurs only when the uj moves with respect to 
the medium. 

The expression for the decay width reads 



r = ! 3 ivi f d3p 

2 n s ^ 2q J (2tt) 3 2 Po 
x\M\ 2 (2TT) 4 S^{q- P - P r ), 



d 3 p' 



(2tt)32^ 
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where the factor | is the symmetry factor when the decay 
proceeds into two neutral pions, the factor of 3 accounts 
for the isospin degeneracy of the final pion states (i.e. 
neutral and charged pions), n s is the number of spin states 
of the uj meson, and denotes the sum over these spin 
states (q , p and p' = q — p are the four-momenta of the 
uj meson, and the two pions, respectively). 

We take the effort to analyze separately the longitudi- 
nally and transversely polarized uj, since the presence of 
the medium results in different behavior of these states. 
Transversely polarized uj has two helicity states (n s — 2), 
with projections s — ±1 on the direction of q, and the 
longitudinally polarized uj has one helicity state (n s = 1), 
with the corresponding projection s = 0. An explicit cal- 
culation yields 



\ <r^ cr* V rpflv 

1^ £ (sf(s) = ~ J - 
s=±l 
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There is another possible structure in the ttNA coupling, 
of the form aj/y' 1 . Our vertex corresponds to the popular choice 
of the off-shell parameter a set to zero. 
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Fig. 2. The in-medium width of the uj meson plotted as a 
function of its 3-momentum |q|. The solid (dashed) lines cor- 
respond to the longitudinal (transverse) mode. The labels (a), 
(a-b) and (a-d) refer to Fig. 1. They indicate the diagrams in- 
cluded in the calculation: (a) u) — a mixing mechanism, (a-b) 
uj — a mixing together with the "triangle" nucleon diagrams, 
(a-d) the full result, including the contribution from the A 
isobar. 



The tensors T^ v and L^ v are denned in such a way that 
they are projection operators. In addition, T^ v q v = 
and L^ u q u = 0, which reflects current conservation, and 
T^ v u v — 0. From relations (||) and (0) in Eq. (g) we find 
that 



\Mi 



\A\ 2 Pll {-T^)p v , 



(8) 



\M L \ 2 = (A*p fJi + B*u„)(-L^)(Ap u + Bu„). 



Note that the value of the coefficient C is irrelevant for our 
calculation. For the diagram (a) we find A = 0, B ^ 0, 
hence this diagram contributes only to the width of the 
longitudinal mode. Diagrams (b-d) have 4 ^ 0, B ^ 0, 
and contribute to the width of both the longitudinal and 
transverse modes. 

As we have said, we evaluate the amplitude M. to lead- 
ing order in the baryon density. This leads to a simplifi- 
cation. The integrals of the form J Q F k 2 dk f(k) arising 

in our calculation are replaced by /(0) J Q F k 2 dk, which 
is proportional to baryon density, ps- Consequently, the 
widths r L ' T ~ p%. 



In Fig. 2 we present our numerical results at the nu- 
clear saturation density, pb = 0.17fm~ 3 . We show F L 
(solid lines) and F T (dashed lines) plotted as functions of 
|q|. The upper part of the plot is for to* = m u = 780MeV, 
i.e. the value of the u> mass is not modified by the medium. 
The lower part is for m* = 0.7™^. In both cases we re- 
duce the value of the in-medium nucleon mass to 70 % of 
its vacuum value, M* = 0.7M, which is a typical number 
at the nuclear saturation density. We also reduce by the 
same factor the mass of the A, i.e. M* A — Q.7Ma, since it 
is expected to behave similarly to the nucleon. The labels 
indicate which diagrams of Fig. 1 have been included. The 
complete result corresponds to the case (a-d). The case 
(a) reproduces the result of Ref. We note that the 
inclusion of subsequent processes of Fig. 1 substantially 
increases the result. All the curves start as q 2 at low |q|. 
The longitudinal width reaches a maximum at |q| ~ a few 
hundred MeV, and the value at the peak is large: 250MeV 
for m* = and lOOMeV for m* = 0.7m w . The trans- 
verse width is strictly zero with diagram (a), less than 
lMeV with diagrams (a-b), reach a few MeV when the 
diagrams with the A are included. Qualitatively similar 
results follow for other choices of parameters. One should 
bare in mind that the effect is proportional to p B , hence 
may be much larger at higher densities. 

Our main conclusions are: 1) nuclear matter induces 
the uj — > 7T7T transitions with large partial widths, 2) the 
widths depend strongly on the three momentum of the u> 
with respect to the medium, |q|, 3) the longitudinal mode 
is much wider than the transverse mode. 

The results obtained mean that in a hadron gas, such 
as created in a heavy-ion collision, the propagation of lon- 
gitudinally polarized u meson is inhibited when the mo- 
mentum |q| is nonzero. This will cause a depletion in the 
population of the uj mesons. Such effects should be in- 
cluded in Monte-Carlo simulations of heavy-ion collisions. 

The important question is to what extent can the dis- 
cussed process influence shape of the dilepton-production 
cross-sections in relativistic heavy-ion collisions. Our re- 
sults can be used to calculate the cross section for the 
7T7T annihilation into dileptons occurring in the uj chan- 
nel. This mechanism has been analyzed for the first time 
in Ref. |2l[ . The calculation of the annihilation cross sec- 
tion requires the knowledge of the same amplitude that 
has been used in the calculation of the omega width, see 
Fig. 3. In Ref. pl[ only the uj — a mixing term was in- 
cluded in this amplitude (diagram (a) of Fig. 1). In our 
present calculation we take into account additional dia- 
grams shown in Figs. 1 (b) - (d). Moreover, we take into 
consideration differences in the propagation of the trans- 
verse and longitudinal modes, which is important due to 
the large difference observed in the widths. The 7T7t annihi- 
lation cross section corresponding to Fig. 3, averaged over 
the incoming pion momenta at fixed total four-momentum 
q = (q°, q), can be written in the compact form 



8n 2 q°m 4 u 



9q 2 {q 2 /4-ml) \g, 



(ml\ 2 r T 



\G L ,\ 2 r L ) 



(9) 
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Fig. 3. Pion annihilation process in the uj channel. The 
blob corresponds to the amplitude shown in Fig. 1, and the 
dot denotes the vector-meson-dominance conversion factor 
eml/(2g u ). 

where 

|G " i|2 = [q 2 -mi- i(r + r^.i)2]2 + m 2 ( r + r^y 

(10) 

is the modulus of the u propagator squared. The widths 
r T and r L should be calculated in the way described 
above with the only difference that the physical mass of 
uj is now replaced by the invariant dilepton mass y/tjp. 
The quantity Jo denotes the width of the u) at vanishing 
three momentum, which is due to other effects, such as 
uj — ► iririr. With In ~ lOMcV |l2) and our values for 
r T ' L we obtain the cross section from Eq. ([)]), which is 
typically a fraction of a microbarn. This is to be compared 
to 3.5/ib from the decay via the p resonance [Ml. Note 
that large widths r T L in Eq. (JoJ) do not increase a, since 
they also appear in the denominator of Eq. (|o|). In fact, 
there are optimum widths r T,L ~ Jb at which the cross 
section is the largest. A further increase of the widths F T ' L 
decreases the cross section. At the point q 2 = and with 
our numbers from Fig. 2 we find that the contribution of 
the longitudinal modes to Eq. (^) is negligible, while the 
contribution from the transverse modes at |q| = 400MeV 
equals OA/ib for m* = m w (where r T = 1.8MeV), and 
1.4/ib for m* = 0.7m u (where T T = 5MeV). We stress 
that the numbers quoted above are almost entirely due 
to the diagrams with the A. Without the processes (b-d) 
of Fig. 1 the dilepton production via mechanism of Fig. 
3 would be about a factor of 10 smaller. In conclusion, 
the process of Fig. 3 may be significant for the dilepton 
production in heavy-ion collisions. 

If the dilepton-production experiments measured the 
three-momentum q of the dilepton pair coming from a 
vector-meson decay, then they should observe different be- 
havior at different values of |q|. Such measurement would 
be very helpful for a better understanding of meson dy- 
namics in the nuclear medium. 

Our last remark refers to the final state interactions, 
which can be important (^6). The pions emitted in pro- 
cesses of Fig. 1 can interact in the final channel. This will 
result in an appropriate modification the u> — > tttt ampli- 
tude. The full analysis of the final-state interactions re- 
quires a model for the 7T7t scattering amplitude, as well 
as solving a Lipmann-Schwinger equation. This is beyond 
the scope of this paper. Note, however, that the diagram 



(a), which includes the intermediate a state, does in fact 
account for final-state interactions. In this process the pi- 
ons form a resonance in the iS-channel, which enhances 
the amplitude. Similar rescattering can also occur for the 
diagrams (b-d). Thus, the final-state interactions are only 
partially included in our analysis. 

We thank Bengt Friman for numerous valuable comments and 
for the suggestion to include the A in the presented analysis. 
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